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Abstract 

n 

A simple model of language evolution, proposed in is characterized by a pay-off in commu- 
nicative function, and by an error in learning, that measures the accuracy in language acquisition. 
In the mean field approximation, this model exhibits a critical coherence threshold, i.e. a minimal 
accuracy in the learning process is required to maintain linguistic coherence. In this work, we 
analyse in detail the effects of different fitness based dynamics driving linguistic coherence and 
of the network of interactions on the nature of the coherence threshold, by performing numeri- 
cal simulations and theoretical analyses of generalized replicator-mutator dynamics in populations 
with two types of structure: fully connected networks and regular random graphs. We find that 
although the threshold of the replicator-mutator evolutionary model is robust with respect to the 
structure of the network of contacts, the coherence threshold of related fitness driven models may 
be strongly affected by this feature. 

PACS numbers: 64.60. Cn, 89.75-k, 87.23. Ge 
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I. INTRODUCTION 



Statistical physics has become a powerful framework to investigate the collective behavior 
of individuals and is playing an increasingly prominent role in quantitative social sciences 
studies. A case in point is opinion dynamics [2j that aims at describing the emergent social 
behaviour, by considering models with simple rules of opinion formation through which 
'agents' update their internal state, or opinion, through the interactions with other 'agents'. 
The interactions are typically local rules that consist in (a) following the majority or (b) 
random neighbour imitation, two simple mechanisms that have been studied for decades 
as models for the dynamics of Ising spin systems, known in the physics literature as the 
Glauber and Voter models, respectively 0,0]. 

Traditional statistical physics models consider particles (spins, agents) interacting (i) with 
all the other particles as analytical solutions are often possible in this mean-field limit or 
(ii) with a number of neighbours located on the vertices of regular lattices in d dimensions, 
the topologv characteristic of crystalline solids. Recently, however, the field of complex 
networks |5|, |6|] paved the way for a better description of social dynamics, by providing 
adequate models for networks of social interactions that are neither well mixed as in (i) nor 
completely regular as in (ii). Since then, numerous studies have considered the evolution of 
opinion models on complex networks, and investigated the effects of the network topology on 
the model's dynamical behavior. In particular, novel, non-trivial behaviour has been found 
for the ordering dynamics of the zero temperature Glauber and Voter models on complex 



networks 
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Language competition may be viewed as a particular case of consensus problems and as 



such has motivated related studies 
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141 ] . Other aspects of language dynamics include 



language change and evolution and language learning. In this context we follow the pioneer- 



ing work of [1] and consider an evolutionary game model where the errors in learning are 
assumed to be the major determinant for language change. This class of models for language 
change are based on the assumption that languages evolve like individuals in a population: 
the fittest survive and spread, the less fit are eliminated. The two driving forces evolution, 
selection and mutation (i.e. language transmission with a bias that favours the fittest or 
the dominant language and errors in the transmission process), are incorporated into the 
replicator-mutator dynamics equations and the time scale for change is generational. 
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In the framework of evolutionary models for language dynamics, the question that arises 
is: how accurately children have to learn the language of their parents in order for the pop- 
ation to maintain a coherent language? The question was answered in a series of papers 



ulatio: 



16] that show that in the strong selection limit a critical threshold, largely deter- 



mined by the error rate of language acquisition, exists for infinite (well mixed) populations, 
irrespective of the number of languages in competition [l| . 

More recently, language games such as the naming game used to model the emergence 



17, Ll8fl have attracted the attention of the 



of language understood consensual lexicon 

physics community H]- This class of models focuses on the horizontal transmission 
and 'creation' of language as a result of peer-to-peer interaction, in contrast with vertical 
;ransmission, the basic scheme of language change in models inspired by biological evolution 
l|. Apart from the time scale for change, which is no longer generational, here the question 
is to establish when the dynamics of a set of interacting agents that can choose among 
several options leads to consensus, or alternatively, when a state with several coexisting 
options, or language diversity, prevails. The discovery of critical thresholds in the naming 
game is arguably the most important contribution of statistical physics in a problem of 
computational linguistics. 

Another class of evolutionary models of languages with analogies with the theories of 



population genetics was proposed by Baxter et al [2l|. The model was solved in the limit 
of a single speaker as well as for multiple speakers in the mean-field approximation and (in 
these limits) it was shown to be related to the model of Abrams and Strogatz |22fl for the 
extinction of languages. 

In this work we follow the view on evolutionary models of language dynamics proposed 
by [l| and consider the investigation of noise induced thresholds for linguistic coherence. 
In other words, we focus on the study (both for deterministic and stochastic versions) of 
the effectivness of the rate of learning errors, in precluding the emergence of linguistic 
consensus. In this framework the coherence threshold is the error rate of language acquisition 
above which a multi-lingual community is stable and below which there is a single dominant 
language. 

In the following, we extend the work of |l| by analysing (i) a family of fitness driven 
models based on the replicator-mutator dynamics, that reduce to the Glauber and Voter 
models in the limit of neutral evolution and (ii) non-trivial networks of interaction. The 
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results detailed below are based on the original Komarova-Nowak [l| model with these 
generalizations, although the time scale for change through learning and selection is no longer 
generational. As in other models of social interaction and opinion dynamics, learning and 
selection occur on a shorter time scale, associated with 'horizontal' (peer-to-peer) interaction. 

The basic assumptions of our family of models are that each individual in the population is 
a speaker of one of two languages 1 or —1, and that an individual may change its language 
through interactions with its neighbours. These interactions follow certain rules, where 
the fitnesses of the individual and of its neighbours determine the probability for language 
change in the absence of errors. In line with the usual replicator-mutator dynamics the state 
update that comes out of these rules is reversed with probability u, that models learning 
errors as the presence of noise in the system coupled to the dynamics. In order to assess 
the robustness of the coherence threshold of the replicator-mutator dynamics we consider 
models with more general fitness driven rules, that reduce to the Voter and Glauber models 
in the limit of neutral evolution and zero noise. The latter are models of spin dynamics, 
used to model the mechanisms of opinion dynamics and cultural evolution, that play a role 
in the evolution of languages on short time scales. 

We find that, in general, in well mixed populations dynamical noise is not effective in 
driving a critical coherence threshold. In models that reduce to the Voter and Glauber 
dynamics in the limit of neutral fitness the noise induced thresholds, separating a dominant 
language regime from the regime where different languages co-exist, become non-critical 
in the mean- field limit. We derive analytical solutions for the coherence thresholds of the 
models in complete and regular random graphs, that are shown to provide a very good 
description of the different types of cooperative behaviour of this family of models. In 
particular, the increase in robustness of the coherent linguistic regime as the number of 
neighbours increases is described quantitatively by the analytical solutions, for all dynamical 
models. 

Finally, we put our results in a more general context and provide a complete classification 
of the threshold behaviour of a family of fitness driven models models that includes a flipping 
rate, or noise uncoupled to the dynamics, instead of the dynamical noise of the replicator- 
mutator equations. 
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II. FITNESS DRIVEN MODELS AND DYNAMICS 

We consider the simplest case of the model introduced in jl| characterized by strong 
selection and two equally fit languages, with no affinity between them. Following [1], we 
consider a population of N individuals, where each individual i speaks one of two languages 
crj = ±1 and define the fitness j\ of i as the number of its neighbors that speak the same 
language, 

i 

where J2'j is a sum over the neighbors of i and 8^1 is 1 if k = I and otherwise. The 
evolution of the language follows two general rules: (i) the language of the fittest individuals 
at a given time step (generation) has a higher probability of being learnt by the population 
in the next time step (selection); (ii) in the process of learning there is a probability of error, 
i.e. a probability that the new generation learns a language with a lower fitness (mutation). 

It was shown in [l| that, in well mixed infinite populations, this model exhibits a critical 
coherence threshold, determined by the rate of learning errors, below which a dominant 
language is established and maintained in the population. In what follows we analyse the 
robustness of the linguistic coherence threshold when other mechanisms of evolution and 
networks of interaction are considered. 

In particular, we consider two generalizations of the replicator-mutator model by in- 
troducing: (i) more general fitness driven dynamics (including additional imitation/social 
pressure mechanisms) and (ii) populations with non-trivial interaction networks. 

We define the social fitness of an individual speaker, F±(i,t), as the total fitness of the 
neighbors of % that speak language ±1 at generation t: 

^±M = £/;(0<W) )± i (2) 
i 

and denote by u the probability of learning errors (0 < u < 1). Unless otherwise stated 
the population has a fixed number N of individuals. The language (of the population) at 
a given time step, t, is characterized by the array {eri(t), (Tzif) ■ ■ - &N(t)}. In the next time 
step the language is determined by the probability that each individual changes its language, 
through the combined effect of the dynamics and learning errors. We have considered three 
fitness driven dynamical models, with learning errors, as detailed below. 
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A. Replicator-mutator dynamics 



In this model the noise, or rate of learning errors, is incorporated in the probability of 
language change. If Gi(t — 1) = 1 the probability for language change (i.e., the probability 
that the outcome of the update rule is <Ji{t) = —1) is given by jl], 



(l-u)F_(i,t-l)+uF + (i,t-l) 
M F„(i,t-l) + F+(i,t-l) ' 1 ' 

while if &i(t — 1) = — 1 the probability for language change (i.e., the probability that the 
outcome of the update rule is <Ti(t) = 1) is given by, 

(l-u)F + (i,t-l)+uF_(i,t-l) 

F_(i,t-l) + F + (i,t-l) " 1 ) 

B. Fitness driven Voter dynamics 

The update rule for this model is inspired in the simplest opinion/imitation dynamics 
model, the Voter model: a speaker changes language if its fitness is lower than a randomly 
chosen neighbour that speaks a different language. In addition, with probability u the 
outcome of the dynamical rule is reversed. The update rule for a speaker i at time t is 
the following: (i) choose one neighbor, j, at random; (ii) if fj(t — 1) > fi(t — 1) then 
(Ti(t) = (jj(t — 1); (iii) if fj(t — 1) < fi(t — 1) then cjj(t) = <Tj(t — 1); (iv) reverse the outcome 
of the dynamical rule with probability u. 



C. Fitness driven Glauber dynamics 



The update rule for this model is inspired in the Glauber dynamics of spin systems, that 
mimicks the effect of social pressure in opinion dynamics. In this model each individual 
tends to adopt the fittest (dominant) language in its neighborhood. The update rule for 
a speaker i at time t is the following: (i) if F + (i,t — 1) > F^(i,t — 1) then <Ji{t) = 1; 
(ii) if F+(i,t- 1) < F-(i,t- 1) then er<(t) = -1; (iii) if F+(i,t - 1) = F_(i,t - 1) then 
<7i(t) = <7i(t — 1); (iv) reverse the outcome of the dynamical rule with probability u. 
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III. SIMULATIONS FOR WELL MIXED POPULATIONS AND FOR REGULAR 
RANDOM GRAPHS 



The models described above were first simulated on complete graphs or fully connected 
networks, i.e., where all individuals are neighbors of each other. On these networks the 
fitness of individuals speaking the same language is identical, in each time step (generation): 
the fitness of individuals speaking +1(— 1) is N% — 1 (N — iVj — 1), where N\ is the total 
number of speakers of 1. 

We start the simulations from a fully ordered system, i.e., at t — all individuals speak 
language +1 (say). The language of the next generation is determined according to the rules 
described above for each model. The language of the population evolves through a large 
number of generations (5000) and the mean value of x = ^ is calculated for each value of 
error rate in learning, u. 

In figure 1 we plot the results of simulations of the three dynamics, for different population 
sizes, N = 100, N = 1000 and N = 10000. Although finite size effects are visible for the 
smaller systems they are negligible for populations of thousands. While we find a critical 
coherence threshold, at u = 1/4, reproducing the results of [l|] for the Komorova-Nowak 
model, the results for the Voter and Glauber fitness driven models are quite different: the 
threshold is shifted to the value of the noise that completely overrides the dynamics, u = 1/2, 
and the fraction of speakers of the dominant language approaches x = 1/2 linearly, at the 
threshold, revealing its non critical nature. 

In order to investigate the effects of the network of interactions on the linguistic coherence 
threshold we simulated the same models on regular random graphs (RRG), where analytical 
results may also be obtained. In RRG networks N nodes are linked at random to a fixed 
number of neighbors, k, without double and self links. The models were simulated on two of 
these networks, for small (k = 4) and large (k = 20) degree. The simulations, for N = 10 3 
and N = 10 4 , start (as before) in the ordered state where all individuals speak +1. For each 
value of u, the language evolves through 10000 generations at the end of which the average 
fraction of speakers of the dominant language is computed. The results for the three models 
are plotted in figures 2,3 and 4. 

Note that the transition on RRG exibhits a critical threshold for all models. The value 
of u at threshold, Uth, increases with the number of neighbours and approaches the MF 
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values (| for the replicator mutator and | for the voter and Glauber dynamics) as the 
number of neighbours tends to infinity. Above threshold, u > u t h, the equilibrium value of 
x corresponds to the coexistence of the two languages, x = \. 

IV. ANALYSIS OF THE MEAN FIELD EQUATIONS 

In order to shed light on these results we proceed to calculate the equilibrium values of 
Ni analytically, in the infinite population limit. Let x be the fraction of speakers of language 
1 (x = Ni/N). The evolution of x is given by, 

x = -xPi^i + (1 - x)P_i-i (5) 

where Pi_>_i and P_i_»i are the rates of change of the two competing languages. 

In well mixed populations Pi_,._i and P-1-+1 depend only on x and can be computed 
exactly for the three models. 

On RRG networks these probabilities are calculated using the following (mean-field) 
assumptions: (i) each of the k neighbours of any site are linked to (k — 1) second neighbors, 
with no loops (uncorrelated links); (ii) the probability that the language spoken at a given 
site is +1 is the average density of speakers of that language, x (uncorrelated densities). 
Within this mean-field approximation, the fitness of each node is a random variable that 
results from the sum of k independent and identical binomial variables. The calculation of 
the transition probabilities Pi-»_i and P_i^i depends on the specific dynamics and proceeds 
in a straightforward fashion. 

Given the symmetry of the models the probabilities © may be written as, 

Pi— i = {l-u)Q(x)+u(l-Q{x)) (6) 
P-i->i = (1 - u)Q{\ -x) + u(l - Q(l - x)), (7) 

where Q(x) is a function that depends on the network of contacts and on the dynamics. 
Substituting ( JoTTl) . (jHJ) becomes 

x = (1 - 2u) (-xQ(x) + (1 - x)Q(l - x)) + u(l - 2x). (8) 

In what follows we discuss the meaning of Q(x) for each dynamical model and calculate 
it for each network of contacts. 



S 



A. Replicator-mutator dynamics 



We find by inspection of and (IS|7j) that for the replicator-mutator model, Q(x) 

(Q(l —x)) is the normalized value of F- (F+) in the neighbourhood of a speaker of language 
1 (-1), 

QW = j 1 ^- 0) 

In well mixed populations, F_ and F+ take their mean values, and Q(x) becomes, 



Substituting (Till into flSJ) yields the evolution equation 

(-x(l - x) + u) (1 - 2x) 
X 2 + (1 — X 

and the (stable) fixed points, x*, are solutions of x = 0. It is straightforward to show that 



X 



| if u > I 



(12) 

\ (l ± v / T Tr 4^) if « < I 
confirming that = ^ is the threshold for linguistic coherence. Furthermore, this threshold 
is critical since the derivative of x* w.r.t. u diverges there. The function (112j) is plotted, for 
x* > |, in figure 1, and excellent agreement is found between the analytical solution and 
the simulation results for large systems. 

The transition probabilities of the replicator-mutator model on RRG are calculated by 
determining the average value of the total fitnesses F + (i) and F-(i) in the neighborhood of 
a given node. Let us consider a node i that speaks +1 with n neighbours that speak also 
+1. The average fitness of one of these neighbours is 1 + {k — l)x and that of the neighbours 
speaking —1 is (k — 1)(1 — x). The average values of F + (i) and F_(i) are then, 

F_ = (fc-n)(A;-l)(l-a:), (13) 
F+ = n(l + (k-l)x). (14) 

The number of neighbours of % speaking the same language, n, is a random variable that 
results from the sum of k random binomial variables, each one taking the value 1 with 
probability x and with probability (1 — x). The function Q(x) in © is then, 

Q(x) = £ B(k, n)x n {l - x) k ~ n -^—r, (15) 



where F_ and F + are given by (I13|I14I) and B(i,j) is the binomial coefficient, 

= (16) 

The evolution equation is obtained by substituting (IT51) in (|Sj). The stable fixed points as a 
function of the noise parameter are ploted in figure 2 for k = 4 and k = 20, respectively. For 
regular random graphs with k = 20 the agreement between the simulation and the analytic 
results is almost quantitative, for populations of the order of a few thousand. 

B. Fitness driven Voter dynamics 

In the fitness driven Voter dynamics an individual that speaks language 1 changes to 
language —1: with probability (1 — it) if the neighbour chosen at random speaks —1 and 
has a higher fitness; with probability u if the neighbour chosen at random speaks 1 or has a 
lower fitness. Thus, Q(x) (Q(l — x)) is the probability to find a neighbour that speaks —1 
(1) and has a higher fitness. In the mean field approximation, Q(x) is 

Q(x) = (l-x)H(x), (17) 

the product of the probability (1— x) of finding a neighbor that speaks —1 and the probability 
H(x) that a speaker of —1 has higher fitness than a speaker of 1. 

In well mixed populations, the probability that a speaker —1 has a higher fitness is 1 (0) 
for x < \ [x > |) implying that, 

Q(x) = (l-x)6(l-2x), (18) 

where Q(z) is the step function: Q(z) — 1 if z > and Q(z) = if z < 0. The dynamical 
equation is obtained by substituting (fT8"|) into © and has stable fixed point solutions, x = 0, 
given by 

'\ \iu>\ 



(19) 

\ (l ± (2a - VI + 4a 2 )) if u < \ 



x 



where a = jz^- Again, the rate of learning errors u th = | defines two regimes: for u < u th 
a dominant language is established and maintained while for u > u th there is coexistence of 
the two equally probable languages. Note, however, that u th = 1/2 is a trivial threshold in 
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the sense that for this level of noise the evolution is totally random, while for higher levels of 
noise the evolution rules actually hinder linguistic coherence. This trivial threshold is non 
critical since, 

dx* 1 , , 

lim — = ±-, 20 

„_»i- du 2' K ' 

the derivative at threshold is finite. The function (TIT?]) is plotted, for x* > |, in figure 1, 
and apart from the finite size effects mentioned previously, quantitative agreement is found 
between the analytical solution and the simulation results. 

To calculate H(x) for RRG, let us consider a node i with a, = +1, and a neighbor j 
with o.j = —1. Using the definition of fitness and the rules of the Voter dynamics we can 
compute H(x) as the probability that j has a number m of neighbours speaking —1 that is 
larger than the number n of neighbours of i speaking +1, 

k-2 



H(x) = J2 B ( k ~ !> n)x n (l - x) 



k—l—n 



X 



n=0 
fc-1 



x B{k-l,m)(l-x) m x k - 1 ~ m . (21) 



m=n+l 



The stable fixed points are calculated using ([H]) with Q(x) given by (11 7112 ip and solving for 
x = 0. In general (x ^ |) they are more easily written in terms of u(x*), 



u[x 



Yi£) jf x * > 1 

2x*-l+V(x*) 2 



n 1 -*') if x * < 1 



(22) 



where V"(x) = s(l — a;)(if(l — x) — H(x)). u(x*) given by fl2^|) is plotted in figure 3 for 
x > 2- The calculated x*(m) is in line w: 
agreement becomes nearly quantitative. 



x > 2- The calculated x*(m) is in line with the results of the simulation and for k = 20 the 



C. Fitness driven Glauber dynamics 

In the fitness driven Glauber dynamics, Q(x) is the probability that a speaker i of +1 
has total fitness satisfying -F_(i) > F + (i). 

In well mixed populations, F_ > F + iff a; < ~ and therefore, 

Q(x) = 6(1 -2x). (23) 
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Substituting in (|SJ) and solving for the stable fixed points we obtain 

x* = I (24) 
\±{\-u) if u < 1 

The rate of learning errors = | is again a trivial threshold that separates two regimes 
as in the Voter driven model. Also as in the Voter model the threshold is non-critical since, 
the derivative of u at threshold is finite and the fraction of speakers of the dominant language 
approaches the threshold linearly. The function (|2"4T) is plotted, for x* > ~, in figure 1, and 
excellent agreement is found between the analytical solution and the simulation results, for 
large systems. 

In order to calculate Q(x) for RRG, we consider in the neighborhood of a given node i 
that is a +1 speaker: (i) n nearest neighbors that speak +1; (ii) ri\ next nearest neighbours 
that speak +1 and share with i a nearest neighbour that speaks +1; (iii) n 2 next nearest 
neighbours that speak —1 and share with i a nearest neighbour that speaks —1. Then the 
total fitnesses are simply given by F_(i) = n 2 and F + (i) — n+rii and Q(x) is the probability 
that n 2 > n + rii, 

k n(k-l) 

Q{x) = ^5(A;,n)a; n (l - x) k ~ n ^ B(n(k - l)^^^ - x) n(k ^ ni x 

n=0 m=0 
(fc-n)(fc-l) 

x £ B({k-n)(k-l),n 2 )(l-x) n2 x {k ~ n){k - 1) - n2 (25) 

The fixed points are calculated using (JHJ) and solving for x = 0. The function u(x*) obtained 
is of the form ([22]) with V(x) = -xQ(x) + (1 - x)Q{\ - x) and Q(x) given by (J25]). u(x*) 
is plotted in figure 4 for x > \ and is found to be in line with the results of the simulations. 
For k = 20 the agreement is almost quantitative. 

Again, the transition exibhits a critical threshold and the value of u t h also increases with 
the number of nearest neighbours on the network approaching the MF value, uth = 1/2, as 
this number approaches infinity. Above threshold, u > Uth, the stable fixed point corresponds 
to the coexistence of the two languages, x = \. 
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D. Coherence thresholds for social fitness driven evolution 



For all the models considered above, the mean field description of the dynamics is given 
by equation (jSJ) which is of the form 

x = (1 -2u)g(x) + u(l -2x) (26) 

with g(l/2) = and g(x) = —g(l — x). For ( 1261) to describe the mean field dynamics of 
an evolution process that selects for the dominant variant of two languages or species with 
the same intrinsic fitness, the additional assumptions are that g(0) = and g(x) < for 
< x < 1/2, so that the coherent states x = and x — 1 are the only stable solutions in 
the absence of noise. The threshold behaviour of this type of models is easily understood if 
we consider the related family 

x = g(x) + u(l - 2x) (27) 

where g(x) has the same symmetry properties as g{x) and noise and dynamics are uncoupled, 
so that u represents a constant rate of random flipping independent of the evolution rules. 
Indeed, ( 12"B1) can be brought to the form (T27j) with g(x) = g(x)(l — 2x)/(l — 2x — 2g(x)) 
through a smooth rescaling of time, provided that g(x) is smooth. From equation (1271) . the 
curve u(x) that relates the rate of random flipping u with the corresponding equilibrium 
density x is 

»(x) = gL. (28) 
If we assume that g(x) is smooth, then given the symmetry 

g{x) = g'(l/2)(x - 1/2) + 0((x - 1/2) 3 ) (29) 

and therefore 



th 



du th 



u~ = u (l/2) = S'(l/2)/2, — =u'(l/2)=0. (30) 

This means that models ( 1271) . ( 1261) always exhibit critical coherence thresholds when g, g 
are smooth. The values of the critical thresholds found in this section are particular cases 
of equation fl30p . which for model (1261) reads 



.th 



u 



1 g'{l/2) du 



tit 



" (1/2) = ITTtm' T x ="'(V2) = o. (3i) 

Equation (1281) also shows that whenever g(x) is a step function with a discontinuity at 
x = 1/2, then model (|2T|) has no coherence threshold: a dominant language persists for 
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arbitrarily large levels of noise, because the right and left limits w(l/2 + ) and u(l/2~) are 
both infinite. The behaviour of model when g(x) is a step function with a discontinuity 
at x = 1/2, as for the fitness driven Voter and Glauber dynamics on the complete graph, 
may be obtained directly from the analogue of equation (l28i) for model ( T26l) 



Then 



«M = 2.-1+ W < 32 > 



w t/l = m(1/2 + ) = u(l/2~) = 1/2 (33) 



independently of g, and 

I* -"' (1/2+ ^ = "55(pv7- (34) 

which is always bounded away from zero. Therefore, this class of models will exhibit a non 
critical trivial threshold at u — 1/2. 

To summarize, the family of models described at the mean field level by equations 
( 1261) . ( 1271) exhibits three types of threshold behaviour. Models with non-smooth density 
dependence transition rates g(x) and dynamically coupled noise exhibit a trivial non-critical 
threshold at the value of u = 1/2 for which noise completely overrides the dynamics. Models 
with non-smooth density dependence transition rates g(x) and dynamically uncoupled noise 
do not exhibit a noise induced threshold, i.e. there is always a dominant language irrespec- 
tive of the level of noise. Finally, in the generic case of models with smooth transition rates 
g(x), g(x) there is a critical threshold, below which language coherence may be established 
and maintained. 



V. CONCLUSIONS 



We have considered different models for the evolution in the presence of noise of two 
languages with the same intrinsic fitness that compete through the selective advantage of 
the language that is perceived by each individual as the dominant language. The language 
spoken by each speaker has for that speaker a social fitness given by the number of its neigh- 
bours that share the same language, and the dynamics driven by evolution rules based on 
this fitness measure will depend also on the interaction network of the population. Starting 
from a state where all individuals speak the same language, mutations or transmission errors 
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act as noise terms that favour the balance of the number of speakers of each language, while 
selection according to social fitness drives linguistic coherence. The coherence threshold is 
the level of noise or mutation rate above which the system evolves to a state where both 
languages are equally frequent. 

From simulations of these models on fully connected networks and on regular random 
graphs, we found that the critical threshold for the replicator-mutator model [l| is robust 
with respect to the network structure, but that the coherence thresholds of related models 
are strongly affected by this feature. In particular, we have found that models with social 
fitness driven dynamics inspired by the Voter and Glauber models, two of the simplest 
models for spin dynamics used in opinion dynamics and cultural evolution studies, exhibit 
different linguistic coherence threshold behaviour, depending on the network of interactions. 
On a regular random graph, these models have a critical coherence threshold, while on the 
fully connected network a dominant language persists up to the level of noise for which the 
evolution rules are totally random. 

We have obtained analytical mean-field solutions for the coherence thresholds on the fully 
connected network and on regular random networks that are in agreement with the results of 
the simulations for the three models, providing a quantitative description of the behaviour of 
the different microscopic rules. We have shown that the noise threshold behaviours of these 
models and, more generally, of evolution processes that select for the dominant variant of 
two languages or species with the same intrinsic fitness, can be understood as well in terms 
of a simple mean field analysis. 
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Figure 
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FIG. 1: Symbols: mean fraction of speakers of language 1 (x) during 5000 generations for different 
values of u and a population N = 10 4 (simulation results in a fully connected network, for an 
initial condition x = 1). Circles: Replicator-mutator dynamics; Triangles: fitness driven Glauber 
dynamics; squares: fitness driven voter dynamics. Lines: fixed points x* from (|T2l [T9l for 
x>\. 
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FIG. 2: Symbols: mean fraction of speakers of language 1 as a function of u for the Replicator- 
mutator dynamics from simulations (10000 generations for an initial condition x = 1.) in a random 
regular lattice with degree k; full symbols: k = 4; open symbols: k = 20; circles: N = 10 3 ; squares: 
N = 10 4 . Lines: fixed points x* from ([8| I15p for x > ^; full line: k = 20; dot-dashed line: k = 4. 
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FIG. 3: The same as in figure [2] but for the fitness driven voter dynamics. The fixed points x* 
were calculated using ([HI 121)170 for x > \ . 
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FIG. 4: The same as in figure [2] but for the Glauber driven voter dynamics. The fixed points x* 
were calculated using ([8l [25]) for x > \ . 
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